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ABSTRACT 


A temperature dependent hard sphere diameter is found from viscosity data for 
sixteen liquid metals. This diameter is used to calculate self-diffusion coefficients using 
the corrected Enskog equation and Stokes- Einstein expression. While both these models 
underpredict the value of the molecular diffusivity, the accuracy of prediction is 
comparable to that obtained using the Protopapas-Parlee hard sphere diameter that is 
obtained from the diffusivity data itself. This scheme is proposed herein has a distinct 
advantage because viscosity data is available for many more systems than the diffusion 
data. The universal applicability of corresponding states model has also been tested using 
various potential parameters; However, the findings are somewhat inconclusive at least at 
this stage. 
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Chapter 1 


INTRODUCTION 

1.1 Introduction 


The study of transport coefficients of liquid metals is important for various scientific 
and engineering reasons. A correct theoretical interpretation of their values would allow 
us to draw conclusions about the nature of ionic and electronic motions and interparticle 
interactions in such liquids. Also, a method of estimating their values would be of great 
use in the design of chemical, metallurgical and other industrial processes involving mass 
transfer by diffusion. 

The process of diffusion is defined as the transport of mass from one region to 
another on an atomic scale. Diffusivity values in liquid state are much higher than those 
in the solid state. In the case of molten metals, diffusivities in the two states differ by a 
factor of 100 to 1000. The high atomic mobility of most metals just above their melting 
temperatures, wherein diffusivities of the order of 10’ 9 m 2 /s, is one of the most 
characteristic properties of liquids. 

The knowledge of self and binary diffusion coefficients in liquid metals is important 
for rational prediction of mass transfer rates in a number of existing and potential 
applications including the pyrometallurgical reprocessing of spent fuels from nuclear 
reactions, and in the production of metallic alloys in liquid phase. Several other processes 
such as corrosion, phase separation, crystal growth, etc., are governed by the rates of 
mass transfer, which, in turn, are primarily dependent on the value of the relevant 
diffusion coefficient. Experimental values of diffusion coefficients are used whenever 
available. Such data are, however, limited because of the difficulties inherent in the 
experimental determination of diffusivity especially at elevated temperatures. 
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Furthermore, it is unlikely that experimental values will ever become available for all 
systems of practical interest and/or at all values of temperature. It is thus often necessary 
to either predict the missing values or to extrapolate the existing values of diffusion 
coefficients; hence the need for a sound theoretical framework for predicting self- 
diffusion coefficients in molten metals is obvious. More reliable and accurate 
experimental data for self-diffusivities of various liquid metals are therefore needed so as 
to test and develop theories. 

Different techniques have been developed for the measurement of diffusion 
coefficients of liquid metals system but theoretical studies have so far been limited to 
those involving a few liquid model approaches. Three distinct theoretical approaches 
have been used in the past to predict the values of self-diffusion coefficients of liquids. 
They are 1) the study of particular models for the liquid state combined with proposed 
mechanism for self diffusion, 2) the application of corresponding states principles, and 3) 
the use of molecular dynamics calculations. 

The most familiar and notable example of the first approach is the Stokes-Einstein 
relationship between the self-diffusion and viscosity coefficients of a fluid. It is derived 
by assuming that the molecular fluid surrounding a diffusing particle can be described as 
a continuous fluid satisfying the usual hydrodynamics equations. Another example is the 
jump diffusion model; in this approach, the mechanism of motion is the jump of single 
atom or molecule from one equilibrium position over an activation energy barrier into a 
neighboring vacant site (hole) in the liquid. 

The second basic approach for the liquid state is hard sphere theory. Hard sphere 
models have long been used to estimate self-diffusion coefficients of liquids in general 
and liquid metals in particular (Ascarelli and Paskin, 1967; Vadovic and Colver, 1970). 
The main reason for this is that the hard sphere model is a good approximation at least for 
the prediction of liquid properties. Further, analytical expressions are available for the 
diffusion coefficient (as well as for other transport properties) of hard spheres. These 
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expressions require, apart from the temperature and density, only one parameter, namely 
the effective hard sphere diameter. Often the hard sphere diameter is used as an 
adjustable parameter chosen so as to fit the predictions of the theory to available 
experimental data on diffusion. Clearly, the predictive capability of such an approach is 
very limited. 

Studies involving a priori determination of hard sphere diameters for transport 
property calculations have been relatively few. Ben-Amotz and coworkers (Ben-Amotz 
and Herschbach, 1990; Ben-Amotz and Wills, 1993) proposed a scheme for determining 
hard sphere diameter from liquid compressibility data. Using this diameter in conjunction 
with Enskog theory, translational and rotational diffusion coefficients in n-hexane were 
predicted with reasonable success. However, given that Enskog theory itself is inadequate 
for hard spheres at higher densities, the need for a most sound theoretical basis is 
obvious. The Enskog theory is based on the assumption that an atom undergoes a series 
of hard-core binary collisions with its neighbors and that successive binary collisions are 
uncorrelated with each other. This approach has been used by Dymond and Alder (1966) 
to develop a theory of diffusion in Van der Waals fluid and by Ascarelli and Paskin 
(1968) to develop a theory for liquid metals. Many other examples of this type of 
approach could be cited, such as the fluctuation theory of Swalin (1959), the free volume 
theory of Cohen and Turnbull (1959) and the significant structure theory of Eyring and 
coworkers (Erying et al., 1964; Breitling and Erying, 1972). 

For Lennard-Jones fluids, use of perturbation theory methods to estimate effective 
hare sphere diameters appears adequate for the prediction of diffusion coefficients 
(Ruckenstein and Liu, 1997). In the case of liquid metals, however, a similar approach 
was not quite as successful (Chauhan et al., 2000). This could be attributed to the long- 
range nature and softer core of the effective potentials of liquid metals. A further 
difficulty with the above approach is its reliance on accurate potentials, which may not be 
readily available. 
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Each of these theories makes implicit or explicit assumptions about the structure of 
a liquid and the nature of diffusive motion. The validity of such a theory must be judged 
partly on the basis of the correctness of its physical assumptions. Agreement between 
theoretical and experimental results is not sufficient to validate a theory, especially when 
it contains one or more disposable parameters. In fact, several of the theories mentioned 
above have proven to be fairly successful in the calculation of self-diffusion coefficients, 
despite the fact that they are based on completely different assumptions about the nature 
of transport process in liquids. I 

The necessity for making detailed assumptions about the liquid state is avoided in the 
second basic approach mentioned above, namely the use of the corresponding states 
principles. In this approach, the reduced self-diffusion coefficient for any liquid metal in 
any thermodynamic state is assumed to be a universal function of the reduced 
thermodynamic properties of that state. Pasternak and Olander (1967) have presented 
such a theory for self-diffusion in liquid metals, which builds upon the approach of 
Chapman (1966). Corresponding states methods are very useful for estimating transport 
coefficients, but they give little insight into the nature of liquids. 

The third basic approach to the study of self-diffusion is the use of molecular 
dynamics (MD) calculations. In this approach, the atoms in a fluid are assumed to interact 
with each other by a certain type of inter-atomic potential. The Newton’s equations for a 
finite number of such atoms are solved numerically, and the solution is used to calculate a 
wide variety of measurable properties, including the self-diffusion coefficient. A number 
of models (Verlet, 1967; Rahman 1964; Alder et al., 1970) have been investigated using 
this method. Molecular dynamics calculations provide a wealth of information about 
structure and atomic motions in liquids, but they are prohibitively expensive in terms of 
computing requirements, and hence, cannot be routinely applied to all fluids of interest. 
However, from the available results of such calculations, it is possible to draw several 
general conclusions, which are useful in developing a general theory of self-diffusion. 
These conclusions can be summarized as follows: 
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1) Dense liquids with realistic interatomic potentials have a structure, which is in 
many ways similar to that of hard sphere fluid. 

2) The self-diffusion coefficients of realistic model liquids are not sensitive to the 
hardness or softness of the repulsive part of the potential, and are similar to those of the 
hard sphere fluid under similar conditions. 

3) For the hard sphere fluid, the physical picture behind the Enskog theory is 
remarkably accurate over the entire range of fluid densities. 

4) The self-diffusion coefficient for hard spheres is very nearly the same as the result 
of the Enskog theory. 

Thus molecular dynamics calculations allow circumventing the problem common to 
most theories of liquids, namely the uncertainty about the validity of the basic concepts. 

1.2 Scope of present work 


In view of the above considerations, it seems worthwhile to pursue an approach that 
retains the simplicity afforded by the hard sphere model but one that neither uses the 
diffusion coefficient data nor depends on the availability of potentials. Instead, the hard 
sphere diameter is determined from the viscosity data, which is generally more accurate 
than the diffusion data. This approach was adopted for two reasons. One, viscosity data is 
available for a large number of liquid metals and hence if this procedure is proved 
successful, it would enhance our capability for calculating diffusion coefficients of liquid 
metals. Secondly, in an earlier work, use of the hard sphere diameter prescription of 
Protopapas and Parlee (1974) in the Stokes-Einstein equation led to a satisfactory 
prediction of diffusion coefficients. It would be interesting to see if the connection 
between diffusivity and viscosity, when established indirectly through the hard sphere 
model has any advantages over the Stokes-Einstein model. 

On the other hand corresponding states model has also been tested with various types 
of potential parameters to obtain a universal curve for viscosity and self-diffusion 
coefficients using Helfand and Rice (1960) equations for corresponding states principle. 
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Chapter 2 


LITERATURE REVIEW 


t 

In this chapter, we will review the literature connected with the theoretical and 
empirical calculations of self-diffusion coefficient in liquids in general and in molten 
metals in particular. It is convenient to classify the available models into the following 
four categories: 

1. Hard Sphere models 

2. Hydrodynamic models 

3. Corresponding States Theory 

4. Other models 

Now, the salient features alongwith their relative merits and demerits are recapitulated in 
the ensuing sections. 

2.1 Hard Sphere models 


This model takes into account the finite size of the molecules but neglects attractive 
forces. Furthermore, molecules are thought to be made-up of hard spheres just like 
billiard balls. For hard sphere molecules, there are no forces between the molecules when 
their centers are separated by a distance larger than a, the hard sphere diameter, but the 
forces of repulsion become infinitely large when they touch, at a separation equal to the 
hard sphere diameter cr. 

The main feature of this class of hard sphere models is that all real fluids behave as 
hard sphere fluids in the high temperature limit. At high temperatures, the real fluid 
equation of state approaches the hard sphere equation of state, and it is no longer 
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necessary to use the experimental density values to get good results. The potential 
function of a hard sphere fluid is expressed by the following equations: 

<(>(r) =0 forr>a (2.1a) 

= oo forr< cr (2.1b) 

Different authors have proposed a range of models for the estimation of viscosity and 
diffusion coefficients based on hard sphere theory. However, we will discuss some 
models, which have been found to be most useful in predicting the values of viscosity and 
self-diffusion coefficient of liquid metals. 

2.1.1 An equation based on hard-sphere theory 


An expression for the shear viscosity of a dense fluid of non-attracting hard spheres 
due to Longuet-Higgins and Pople (1956) has been discussed by Faber (1972) and is 
expressed in terms of the packing fraction 17: 


H = 3.8 x 10 


-8 


(MT) 


1 - 


2 

VJ 


(1 -rif 


( 2 . 2 ) 


Similarly Faber (1972) has given an expression for self-diffusion coefficient of liquid 
metals: 


D = 4.9x10 


-10 


(T} 

yMj 


(h ) 3 


(2.3) 


l-l 


where Mis the atomic weight (kg/g-atom), V is the atomic volume (m 3 /g-atom), 77 is the 
packing fraction, ju. is the viscosity of the liquid metals (Pa-s) and D is the self-diffusion 
coefficient of the liquid metals (m 2 / s) . 

Further, Faber (1972) has also calculated these properties at melting point by taking 
packing fraction 77 to be equal to 0.45. Then equations (2.1) and (2.2) lead to: 
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(2.4) 




= 1.6x1 O' 7 



f T ^ 


D m =4.0x10 


-10 


yMj 


V_ 


(2.5) 


where subscript ‘m’ represents melting point conditions. . 

Iida and Gutherie (1988) showed that Equation (2.3) underestimates the value of p m 
by a factor of the order of 0.3. 

2.1.2 Enskog Hard Sphere model 
1. Equation for viscosity 


Based on the hard sphere theory, Enskog (Chapman and Cowling, 1970) derived an 
equation for viscosity. In this formulation, viscosity is expressed in terms of the packing 
fraction and pair distribution function at contact of hard spheres g (cr) as follows: 


V- 


-t- + 3.20 r, + 12.18 V ‘g(o) L 
,£( CT ) 


( 2 . 6 ) 


and a = — ( junkT)* (2.7) 

where m is the atomic mass, k is the Boltzman’s constant, or is the hard sphere diameter 
and tj is the packing fraction. 

2. Enskog Equations for self-diffusivity 

The Enskog theory (Chapman and Cowling, 1970) assumes that an atom undergoes a 
series of hard-core binary collisions with its neighbors, and that successive binary 
collisions are uncorrelated with each other. This approach simply relates the diffusion 
coefficient with pair distribution function, which in turn is related to an appropriate 
equation of state. This simple model has the advantage that exact calculation of transport 
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coefficients at a given temperature and density can be obtained with only one adjustable 
parameter namely, the hard sphere diameter. Though, the correct estimation of the hard 
sphere diameter and hence the pair distribution function of the hard spheres in contact is a 
difficult task which some times leads to the failure of the theory. 

The self-diffusion coefficient can be computed using the Enskog theory as: 


D = _J_ 

D o g{°) 


( 2 . 8 ) 


where D 0 is given as. 


A = 


8n <j ' 


kTy 
\_nrn J 


the packing fraction, rj, is given by the following equation: 

77=— nff 3 
6 


(2.9) 


( 2 . 10 ) 


n is the number density and is given by: n = — 

m 


( 2 . 11 ) 


where p is the atomic mass density. 
Equation (2.8) can now be rewritten as 


D = 


2(Z-1) 


^ m J 


( 2 . 12 ) 


where r is the hard sphere radius = a/2. 

For a hard sphere fluid the radial distribution function, g(a) is related to the 
compressibility factor, Z, by the following equation: 


Z = 1 + 4q g(a) = 


P 

nkT 


(2.13) 


where p is the pressure of the hard sphere fluid. 

The compressibility factor, Z, can be calculated using the equation of state. 

Percus and Yevick (1958) have given an equation for the compressibility factor Z as 
follows: 
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z = 


1 + 77 + f ] 

& - *if 


2 


(2.14) 


Carnahan and Starling (1969) have given an equation of state, which is superior to any 


existing hard sphere equation of state including that of Percus and Yevick (1958), and is 
given by: 


Z = 


1 + 


n ± n^vl 
(1 - nf 


(2.15) 


2.1.3 Vadovic and Colver approach 


Ascarelli and Paskin (1967) calculated the values of self-diffusion coefficient for 


liquid metals by using equation (2.12) with a correction factor 0.73 


r ri ') 

I m 

V J 


. Vadovic and 


Colver (1970) recognized that at atmospheric conditions, the repulsive hard sphere 
pressure must be balanced approximately by the attractive Van der Waals pressure. They 

Z kT 

evaluated the Van der Waals constant (a) at melting point by using a —— - — — . Further, 

tn 

they calculated the values of diffusion coefficients for nine liquid metals at melting point 
using both equation (2.14) and (2.15) with the equation (2.12). In equation (2.12), they 


used a correction factor 0.73 lm / 


V) 


and the hard sphere radius r = 


3 rj m M 
4 XP m n N J 


They 


concluded that the agreement of all calculated values with the melting point experimental 
data was remarkable, particularly as pointed out by Ascarelli and Paskin(1967), the error 
of the data is generally on the order of 10 percent. Further, they concluded that Carnahan 
and Starling equation of state yields the stronger temperature dependence than that 
stipulated by Ascarelli and Paskin (1967). 
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2.1.4 Protopapas and Parlee approach 


Protopapas et al (1973) have presented a semi-empirical approach for the 
calculation of self-diffusion coefficient of liquid metals. The basic assumption is that 
the self-diffusion coefficient of a liquid metal is equal to that of an appropriate hard 
sphere fluid. The hard sphere diameter (Protopapas and Parlee, 1974) is dependent 
upon temperature, and they developed a method for estimating this temperature 
dependence by exploring the relationship between the diameter and the interatomic 
potential energy function of the liquid metal was presented. The scheme for calculating 
the hard sphere diameter developed by Protopapas et al (1973) is given by: 

( T 

cr = 1.126cr ra 1-0.112 — (2.16) 

\T m J 

where a m is the value of a at the melting point and is given by: 



1 



(2.17) 


In the above equation Protopapas et al (1973) have taken r\ m = 0.472. On further 
simplification, Eq. (2.15) becomes as: 



(2.18) 


Protopapas et al (1973) have used Eq. (2.18) in the Enskog equation (2.8) and 
obtained an equation for self diffusion coefficient as: 


D — c C AW ( rj 


tvRT iz 
M , 


(1-7) 3 

877(2-//) 


(2.19) 


where Caw(j ]) is the correction factor. They computed self-diffusion coefficient at 
melting points and compared with the experimental values for a series of liquid metals. 
A good agreement with experimental results with an average difference about 2 % was 
reported. 
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Khanna et al (1993) have calculated the value of the self-diffusion coefficient of 
liquid alkali metals and its temperature dependence using the Enskog hard sphere 
model. In their study, they tested the suitability of the Enskog theory for self-diffusion 
coefficient for various systems of different range of densities. They found that the pair 
distribution function of the hard spheres in contact, g(cr) depends upon the model and 
its parameters. 

2.2 Hydrodynamic Model 

2.2.1 Stokes-Einstein Equation 


Stokes-Einstein equation is derived based on the assumption that the molecular fluid 
surrounding a diffusing particle can be described as a continuous fluid satisfying the 
hydrodynamic equations. Therefore for a single particle moving in a continuum or for a 
single heavy particle moving in a molecular medium of lighter particles. Brownian 
motion theory predicts that the diffusion coefficient is inversely proportional to the 
frictional coefficient (Q 


D = 


kT_ 

% 


( 2 . 20 ) 


where E, is the frictional coefficient. 

Stokes derived an expression for the frictional coefficient E, for the case of a 
spherical particle of radius R, moving in a continuum having a coefficient of shear 
viscosity jj, as: 


£ = 6k^R 


' 2 ^ 
1 + ^ 

m 


i+ 


la 

PR 


( 2 . 21 ) 


where /3 is the ‘coefficient of sliding friction’ between the diffusing particle and its 
surrounding continuum. There are two possible limiting conditions, namely j3 = oo and 


j3=0. 


12 



a) (5 = <x>: 

When the diffusing particle is large, compared with the particle of the medium, j3 
equals to infinity, i.e. there is no slip between the moving sphere and the fluid in contact 
with it. For this no slip condition, the frictional coefficient becomes: 

§ = 6n/jR (2.22) 

b) 15 = 0: 

When the radius of the diffusing particle approximately equal to that of the 
medium then f5 becomes zero and in this case the frictional coefficient becomes: 

% = AnixR . (2.23) 

In the above two extreme cases equation (2.20) may be written as 


D= kT 

when p = oo 

(2.24) 

6%jjR 

d. kT 

when P = 0 

(2.25) 

4n/±R 


From the above two equations, one can show that the quantity should be constant, 

provided that the molecular diameter is independent of temperature. In general, Stokes- 
Einstein equation can be written as: 

kT 

D = where 4 < a < 6. (2.26) 

a nfxR 

Clearly, when considering the phenomenon of self-diffusion, intuitively it appears 
that a = 4 would apply. In fact, sufficient evidence is now available in the literature 
illustrating the applicability of Eq. (2.26) to both self-diffusion and inter-diffusion 
coefficients in molecular non-polar substances provided a judicious choice is made of 
radius R. For a macroscopic solid sphere, there is no ambiguity in the assignment of 
particle radius, R. For a molecule or atom in the case of liquid metals; however, R must 
be interpreted as some average molecular radius. Thus for ordinary liquids, some 

i 
3 

, where n is the number density of the pure substance in 
liquid or solid state at low vapour pressure, while others have used the procedures 


investigators have used R = — 

In 
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outlined by Edwards (1970) to evaluate R. in many cases, all these schemes yield values 
close to each other. 


2.3 Corresponding states model 


According to corresponding states principles, the reduced viscosity or 
diffusivity of similar substances in any thermodynamic state should be a universal 
function of the reduced thermodynamic properties of that state such as reduced 
temperature and pressure if the molecules are simple (monoatomic or spherical) and if the 
potential energy of interaction can be represented by a universal two parameter function 
of the type: 






(2.27) 


K°j 


where e and cr are the energy and distance coordinates of the function minimum. 
Chapman (1966) has obtained an equation for viscosity by using Bom and Green (1947) 
theory which is written as: 


/•* 


An n 2 5 


3 kT 


m 


dr 


[l - g(r)}lr 


(2.28) 


where A is a constant and 8 is the distance parameter. 


Chapman (1966) postulated the reduced viscosity to be a function of the reduced 
temperature T* and reduced volume V, defined as follows: 


M'=f(r,v) 


(2.29) 
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(2.30) 


. juS 2 N 


(MRT)i 


r =— 

8 

(2.31) 

r= 1 
»<5 3 

(2.32) 


The corresponding states principle is further applied on the assumption that g(r) 
should be a universal function of x = r/8 and T*; when this change of variable is executed, 
and when the reduced quantities are substituted in Eq. (2.28), the final form obtained as 
follows: 

ju'CF*) 2 = G(T*) (2.33) 

where G(T*) = C JV [l ~ g(*T )}& 

where C is a constant. 


Thus ju*(F*) 2 is a universal function of T*, and when plotted against 1/T* should 
generate a universal curve for all liquid metals. This is in fact the case. In order to 
establish a corresponding state plot of the function 2 vs. T*, Chapman (1966) 

assumed that all liquid metals obey the same function $>). He used Lennard-Jones 
potential parameters for Sodium and Potassium determined by Ling (1956) and 
experimental viscosity data, then plotted /u*(V*) 2 vs. 1 IT*, using Goldschmidt atomic 
diameter (Smithells, 1962) for 5. Then from this plot Chapman (1966) evaluated e for the 
remaining liquid metals. Furthermore, with the obtained values of e, a relationship 
between s/k and melting point temperature T m was obtained as: 


f = 5.2T m (2.34) 

k. 
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The generalized plot for all liquid metals obtained by Chapman (1966) is shown in fig 2. 



i/r 


Fig 2 Reduced viscosities of liquid metals and their dependence on reduced temperature. 


This universal curve is fitted into the following form: 

K-(r)'- 


0.3501 ( 0.3056^ 

-exp — 

A T J 


(r): 


(2.35) 


Helfand and Rice (1960) have also shown that by assuming a pair potential 




(j>(r) = s(f>* — , where e and 6 are characteristic energy and distance parameters and <j) is 




a universal function of r/8. The reduced viscosity and the reduced diffusivity can be 
defined as: 

ju<5 2 


» = 


(me) 


(2.36) 
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( 237 ) 


s 2 5 

Pasternak and Olander (1967) also correlated the self-diffusion coefficient and viscosity 
data using the corresponding states principles. The reduced viscosity and reduced 
diffusivity coefficient are defined as: 



where R is the universal gas constant. 


(2.38) 


(2.39) 


Pasternak and Olander (1967) used the values of e and Goldschmidt atomic diameter 
5 as obtained by Chapman (1966). He plotted D*/(V*) 1/3 vs. 1/T* for ten liquid metals. 
All the data points were seen to lie on the same line with a spread of ± 20 %. 


2.4 Other models 

2.4.1 Andrade Equation (Semi-theoretical) 

Andrade (1934) regarded the atoms in the liquid state at the melting point as 
executing vibrations about equilibrium positions with random directions and periods just 
as in the solid state. Furthermore, assuming that the transfer of momentum of atomic 
vibrations from one layer to neighboring one produces viscosity of liquid and he gave the 
following equation for viscosity at the melting point. 
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(2.40) 




1.x 10 



This equation reproduces experimental data for liquid metals at their melting points with 
comparatively good accuracy. Furthermore, Andrade (1934) has presented an expression 
for temperature variation of viscosity. 


juo 3 = ^4exp 




(2.41) 


where u is the specific volume and A and c are constants to be determined from 
experimental data. This equation has been found to fit experimental data with high 
degree of accuracy. 


2.4.2 Empirical models 


The temperature dependence of viscosity data in liquid metals can be approximate by an 
Arrhenius type relation and is given as: 


f 

\x = A exp 

v 


Ejl 

RT 


\ 

J 


(2.42) 


where A and are constants. 

Grosse (1961) has obtained the equation for the constant A using Andrade’s formula for 
the melting point viscosity. This constant then becomes: 


A = 5.7 x 10' 5 


(¥T m ) 


VJ exp 


(H., 


(2.43) 


K RT m J 


lida and Gutherie (1988) have reported expressions for H^, which relates rather well 
with the melting temperature T m . 

For normal metals: (h ^ ) b =1.21 T„ 2 (2.44a) 

For semi-metals: {H h = 0.75 Tj 2 (2.44b) 
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Several investigators have presented empirical equations for self-diffusion in liquid 
metals. They have represented self-diffusion coefficient in the form of Arrhenius 
equation. One of such forms has been reported by Iida and Guthrie (1988) and is given 
by: 


D = D o exp 


■H r 


RT 


(2.45) 


where Ho and D 0 are constants. Ho is called the apparent activation energy and it can 
be evaluated by using the following equations: 

For normal metals: (H d ) n = 2.5 Tj 5 (2.46a) 

For semi-metals: (H d ) s = 2.00T^ (2.46b) 


Larsson et al (1972) have reported that self-diffusivity in liquid metals at their 
melting point can be approximated as: 


D 


m 


6 V M 


V 3 J3 


(2.47) 


The factor /3 is expected to be nearly constant for all liquid metals, approximately 
equals to 0.1. Iida and Guthrie (1988) have calculated the values of self-diffusion 
coefficient using this equation and compared with the experimental data. They 
concluded that this empirical equation reproduces the experimental data at melting 
temperature with good agreement. 
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Chapter 3 


ANALYSIS AND VALIDATION 


This chapter deals with the models, which have been used for the calculation of self- 
diffusion coefficients of liquid metals as a function of temperature. 

3.1 Expression for viscosity 


The Enskog hard sphere equation for viscosity is given by: 


Mb = 


— j-r + 3.20 rj + 12.18 7 2 g(cr)U 0 

,#( CT ) J 


where ju o = 


16jcg 2 


(rnnkT ) 


(3.1) 

(3.2) 


Molecular dynamics calculations have shown that Enskog’s expression for a pure 
hard sphere fluid is not exact, since which do not take into account velocity correlations 
in the dense regime. A correction factor, due to Dymond (1976), as given by Castillo and 


Villaverde (1991) for shear viscosity as follows: 

= 1.02 + 10.6l(n -0.495) 3 +247.49(n* -0.813) 3 n> 0.813 (3.3) 

C„ = 1.02 + 1 0.6 l(>T-0.495) 3 0.593 > n> 0.813 (3.4) 

=1.02 0.593 >n (3.5) 

where n* is the reduced density and is given by no*. 

The corrected Enskog equation for viscosity is then, \x corr = (3 .6) 



3.2 Equations for self-diffusion coefficient 


We have calculated self-diffusion coefficient based on the corrected Enskog equation 
and also the Stokes-Einstein hydrodynamic model. In these calculations, we used the 
Protopapas and Parlee (1974) hard sphere diameter (PAP) and the diameters calculated 
from experimental viscosity data. 

3.2.1 Enskog Equation 


The Enskog expression for self-diffusion coefficient is given by: 


^ JD, 

II 

o 

(3.7) 

, 3 ( kTV 

where D o = 

8na \iurn 

(3.8) 

The final form of the Enskog equation is: 

_ 5 ( kry 

Ug — r- 

8nG 2 g(a)\Ttm) 

(3.9) 

where p is the atomic density and is given by: 

P“P.+A(r-r.) 

(3.11) 


where p m and T m are the density and temperature at melting point and A is a 
constant. 

The pair distribution function, g(a) at contact is a function of packing fraction, and 
can be calculated from the compressibility factor, Z. The compressibility factor can be 
calculated from the Carnahan and Starling (1969) and Percus and Yevick (1958) 
equations of states. Since the equation of state due to Carnahan and Starling (CS) (1969) 
is more accurate than that of Percus and Yevick (1958), the CS equation of state will be 
used here. Therefore the pair distribution function at contact can be written g(cr) as: 

( 312 ) 

4rj 
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where Z 


(3.13) 


1 + t] +jf-rf 

(! - vf 

The final form of the pair distribution function at contact, g(a) is: 

Molecular dynamics simulations (Alder et al., 1970) on hard spheres have 
demonstrated the inadequacy of the Enskog results, especially at higher densities. In 
order to facilitate the analytical evaluation of diffusion coefficients, several expressions 
have been developed in which the ‘corrected’ value is expressed as: 

D‘°"=D E f{n) (3.15) 

In this work, we use an expression for f(n ) developed recently by Ruckenstein and 
Liu (1997) and is given by: 


/(«’)= l + 0.94605(n*) 1 5 +1.4022(«*) 3 -5.6898(«*) 5 +2.6626(n*) ? . (3.16) 

3.2.2 Stokes - Einstein equation 


When the radius of diffusing particle is approximately equal to that of the 
surrounding medium, the sliding friction coefficient ]3 becomes zero and in this case the 
Stokes-Einstein equation becomes: 




kT 

AnfiR 


(3-17) 


In the above equation R is the radius of spherical particle moving in a continuum. 

We have calculated self-diffusion coefficient using Protopapas and Parlee diameter 
(PAP) (1974) and with the diameters calculated from the experimental viscosity data. 
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3.3 Effective hard sphere diameter 


The atomic diameter of elements is a parameter of fundamental importance for any 
consideration at the molecular level, dealing in particular with the models of hard 
spheres. The hard sphere diameter of a liquid metal, at a temperature T, understood to be 
the statistical minimum center to center distance of two atoms of metal during a two-body 
collision at temperature, when the metal atoms are assumed to be perfect spheres. 

Many theories, models and approaches have been proposed for the evaluation of 
hard sphere diameters of liquid metals. Dymond and Alder (1966) have shown that the 
hard sphere diameter may be derived as a function of temperature from equilibrium data 
on the basis of the postulate that repulsive and attractive pressure terms in the Van der 
Waals equation of state counter balance each other. Ashcorft and Lekner (1966) have 
derived the temperature dependent hard sphere diameters. In their approach, the 
theoretical Percus-Yevick structure factor is fitted to the X-ray radial distribution 
function of the metal, and the packing fraction q of the best fit yields the hard sphere 
diameter at the temperature of the X- ray experiments. 

A commonly used method is one based on the application of perturbation theory. 
Three distinct prescriptions are due to Barker and Henderson (BH) (1967), Weeks, 
Chandler and Anderson (WCA) (1972) and Lado (1984). All these approaches are based 
on splitting of the actual potential into a harshly repulsive reference potential and a 
weaker long-range perturbation. The three approaches differ in the nature of the splitting 
and also in their prescriptions for the effective hard sphere diameter and have been 
extensively applied to Lennard-Jones fluids with varying degree of success. All these 
three approaches have been discussed somewhere else (Chauhan et al., 2000). 

A semi-empirical approach for derivation of the hard sphere diameters of liquid 
metals has been proposed by Protopapas et al (1973) assuming perfect collisions of atoms 
moving on parabolic potential wells. Their expression includes two constants, the values 



of which were determined from experimental self-diffusivity data. The proposed equation 
for the hard sphere diameter is: 


cr =1.288 


M 


1 

V 


'mj 


f 

1 - 0.112 

V 



(3.19) 


where M is the atomic weight and a is in A 0 units. 


We have calculated the hard sphere diameters from experimental viscosity data. 
In this approach, the corrected Enskog equation for viscosity (3.6) was used. The 
procedure for calculation of hard sphere diameters from experimental viscosity data is 
shown in detail in APPENDIX A. The values of diameter thus obtained are fitted to the 
form: 

a = A + BT + CT 2 (3.20) 

where A, B and C are constants. In the above equation, cr is in A 0 and T is in K. 

Thus with the calculated diameter, we have calculated self-diffusion coefficients 
using the corrected Enskog equation (3.15) and Stokes-Einstein expression (3.17). 


3.4 Corresponding states model 

In this approach, an attempt has been made to test the universal applicability of 
corresponding states equations. Though there are many such approaches available in the 
literature, as discussed in chapter 2, we have chosen the one due to the Helfand and Rice 
(1960). According to their approach, the reduced properties can be expressed as: 



. Dm 1 

s*8 


(3.21) 


(3.22) 



r =— 

s 


(3.23) 



(3.24) 


In the above equations e and 5 are energy and distance parameters respectively. In 
this case, e becomes potential well depth, m is the mass of atom and the superscript 
refers to reduced parameters. 


3.4.1 Effective pair potentials 


In order to calculate the potential depth parameters and atomic distanced, we need 
effective pair potentials for liquid metals. The effective interparticle potential of liquid 
metals is widely modeled as exhibiting a short-range repulsion and long-range 
oscillations. Most of the potentials have been obtained from structure factor data. There 
are, however, several intrinsic difficulties associated with this procedure. First, since the 
structure factor data is relatively insensitive to the details of the potentials, one requires 
very accurate data to obtain an accurate potential. Such data is generally difficult to 
obtain for liquid metals. Furthermore, in much of the early work, these potentials have 
been obtained through route of integral equation theories that are only approximate at 
best and have in fact been since shown to be quite inaccurate. 

Johnson et al (1964) (JHM) have shown that the inversion of experimental 
scattering data by means of the Bom - Green (B-G) (1947) and Percus- Yevick (P-Y) 
(1958) integral equations yields a long-range oscillatory (LRO) pair potentials for a 
number of liquid metals. Their procedure (JHM) was to deduce the effective pair 
potentials from a knowledge of radial distribution function. The JHM long-range 
oscillatory potentials have the form: 

(r 

(p{r) = A — 

\ r 

where A is the amplitude, r 0 is the position of the first node, /3 is the phase, a is an 
attenuation factor and k is given by n/ wave length of oscillations. 


cos(2k{r + /?})exp(- or) 


(3.25) 
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Recent theoretical work (Hausleitner, 1988; Bari et al., 1991; Bari and Jorder, 1992) 
on liquid transition metals has used the effective pair potentials of Wills and Harrison 
(WH) (1983), which are based on a separate treatment of s-p and d states and, in addition, 
take in to account the effect of s-d hybridisation. According to Wills and Harrison (1983), 
the effective pair potential (EPP) is made up of an ^-electron contribution, (jk(r), a d - band 
contribution term, <j^(r) and the repulsive contribution, <j) r (r) from the shift of J-band 
centre, given by: 


<K r ) = <Ps(r) + <k(r) + <t>r(r) 
The individual terms can be given by: 

<t> s ( r ) = Z W cosh (Kr c )— 

r 


<t>r( r )= Z d 


225 h 2 rf 
n 2 m e r s 


(3.26) 


(3.27) 

(3.28) 




fi-Al 

"12 V 

( 28.1^ 

10. 


v K > 




(3.29) 


m e r 


where Z s is the number of s electrons, K is the Thomas-Fermi screening constant and 


is given by K ■■ 


4 K f ^ 


\™o j 


(3.30) 


In the above equation K F is the Fermi wave number and is given by 


K,= 


9xZ s 
V 4r 0 3 ; 


(3.31) 


where r is the interatomic distance, r 0 is the atomic radius, r c is the pseudo potential 
core radius, r d is the Estate radius, Z is the valency, Z d is the number of J-electrons, n is 
the nearest neighbour co-ordination number, a 0 is the Bohr atomic radius, m e is the mass 
of electron and e charge of electron. 

Z d can be calculated by using the following equation: 

Zs=Z-Z d (3.32) 
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Table 3.1 Liquid metals and their Potential parameters. 


Metal 

Potential parameter 

Potential type 

Source 

e(K) 

g(A°) 

Ag 

3325.6 

2.88 

Lennard-Jones (m-n) 

Zhen and Davies (1983) 

Cd 

1561 

3.112 

Lennard-Jones (m-n) 

Zhen and Davies (1983) 

Co 

1108 

2.78 

Wills and Harrison 

Wills and Harrison (1983) 

Cs 


5.2 

Long-range oscillatory 

■HKSBSEEiBESEHHi 


5.19 

Lennard-Jones (6-12) 

Sengupta et al (1977) 

385.5 

4.76 

Pseudo 

Ranganathan and Pathak (1994) 



2.5 

Lennard-Jones (m-n) 

Zhen and Davies (1983) 

Fe 

1892 

2.85 

Wills and Harrison 

Wills and Harrison (1983) 

Hg 

597.85 

5.20 

Long-range oscillatory 

Johnson et al (1964) 

K 

712.7 

4.54 

Lennard-Jones (m-n) 

Zhen and Davies (1983) 

493.37 

4.51 

Long-range oscillatory 

Johnson et al (1964) 

421.4 

4.11 

Pseudo 

Ranganathan and Pathak (1994) 


4.64 


Sengupta et al (1977) 

Li 

713.6 

3.00 

Lennard-Jones (m-n) 

Zhen and Davies (1983) 

447 

3.21 


Johnson et al (1964) 

567.1 

2.73 

Pseudo 

Ranganathan and Pathak (1994) 

371.37 

3.15 

Lennard-Jones (6-12) 

Sengupta et al (1977) 

Na 

798.9 

3.66 


Zhen and Davies (1983) 

516.6 

3.25 

Long-range oscillatory 

Johnson et al (1964) 

421.4 

4.12 

Pseudo 

Ranganathan and Pathak (1994) 

510.633 

3.82 

Lennard-Jones (6-12) 


Ni 

758.12 

2.74 

Wills and Harrison 

Wills and Harrison (1983) 

Pb 

2215.9 

3.48 

Lennard-Jones (m-n) 

Zhen and Davies (1983) 

1126 

3.30 

Long-range oscillatory 

Johnson et al (1964) 

Rb 

610.7 

4.85 

Lennard-Jones (m-n) 

Zhen and Davies (1983) 

377.3 

5.05 

Long-range oscillatory 

Johnson et al (1964) 

402.2 

4.41 

Pseudo 

Ranganathan and Pathak (1994) 

464.21 

4.97 

Lennard-Jones (6-12) 




3.41 

Lennard-Jones (m-n) 

Zhen and Davies (1983) 

Zn 

1681.5 

2.76 

Lennard-Jones (m-n) 

Zhen and Davies (1983) 
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Fig 3.1 Schematic representation of long-range 

oscillatory potentials (Johnson et al., 1964). Fl g 12 Schematic representation of pseudo 

potentials (Ranganathan and Pathak, 1994). 




Fig 3 .3 Schematic representation of WH Fig 3 .4 Schematic representation of Lennard- 

potentials (Wills and Harrison, 1983). Jones potentials (Zhen and Davies, 1983). 
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Many workers have used inter-atomic potentials in the calculation of transport 
properties of liquid alkali metals. Sengupta et al (1977) have calculated self-diffusion 
coefficient of alkali metals using Lennard- Jones (L-J) potentials. They listed energy and 
distance parameters. Ranganathan and Pathak (1994) have used pseudo potentials of the 


form e<f> 



in the calculations of self-diffusion coefficients. 


Zhen and Davies (1983) have developed a new approach' to calculate the Lennard- 
Jones n-m potential energy parameters for metals using experimental cohesive energy 
data. They employed this approach because Lennard-Jones 12-6 potential is not a good 
representation of the potential for metals. One way to improve the agreement between the 
potential model and the experimental data is to use the generalized form of the 
interatomic potential (the Lennard-Jones n-m potentials), 



(3.33) 


where r 0 is the equilibrium nearest distance between the atoms. This generalized 
potential has four adjustable parameters m, n, e and r 0 . 


The schematic representations of these potentials are shown in the following figures 
(3.1- 3.4). The reduced viscosity and reduced diffusivity were calculated using these 
potential parameters and are listed in the Table 3.1. In the calculations of reduced 
viscosity and reduced diffusivity and their volume dependency. Pseudo potentials 
(Ranganathan and Pathak, 1993), Lennard-Jones (m-n) (Zhen and Davies, 1983) 
potentials were used. These results are presented and discussed in the next chapter. 
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Chapter 4 


RESULTS AND DISCUSSION 


4.1 Data accuracy 


Reliable experimental values of both viscosity and self-diffusion coefficients 
encompassing varying temperature ranges were identified from the literature for sixteen 
liquid metals. It is also worthwhile to address the issue of the accuracy of such 
measurements. While experimental uncertainty varies from one technique to another and 
is strongly dependent on the purity of sample. It is generally believed that the available 
diffusion data are accurate to within 5-7% at best and about 20% at worst. Similarly, the 
accuracy of the viscosity data is also dependent on the technique and temperature but 
these are believed to be accurate to within 5% at best and about 10-12% at worst. 

4.2 Hard sphere diameter 

The hard sphere diameter, a, was calculated using from the experimental viscosity 
data via the corrected Enskog expression. The list of sixteen metals, the sources 
containing viscosity data and the temperature range over which the data were used in the 
calculation of hard sphere diameters is shown the Table 4.1. The temperature dependence 
of hard sphere diameter, a, was fitted to the form 

a(T)=A + BT + CT 2 . (4.1) 

The parameters A, B and C were determined for each metal. The best values of 
these parameters along with their correlation coefficients are listed in the Table 4.2. 
Viscosity data for temperatures within +10% of the melting temperature were not 
considered for this purpose. 
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Table 4.1 Liquid metals and their sources for viscosity and diffusivity data 


S.NO 

Metal 

T m 

(K) 

Viscosity 

Temperature 

Range 

(K) 

Diffusivity 

Temperature 

Range 

(K) 

1 . 

Ag 

1233.65 

Gebhart and 
Becker (1951) 


Yang et al 
(1956) 

1260-1400 

Hi 

Bi 

544.1 

Yao and 
Kondic (1952) 

560-680 

Chauhan et al 
(2000) 

540-680 

HI 

Cd 

594 

Miller (1952) 

610-750 

Chauhan et al 
(2000) 

610-750 




Andrade and 
Dobbs (1952) 

325-395 

Chauhan et al 
(2000) 


5. 

Cu 

1356 

Landolt- 
Bomstein 
tables (1969) 

1373-1573 

Henderson and 
Yang (1961) 

1370-1440 

6. 

Ga 

303 

Spells (1936) 

330-365 

Petit and 
Nachtrieb 
(1956) 

330-365 

■EH 

Hg 



253-613 


270-370 




Crawly and 
Thresh (1969) 

440-619 

Lodding (1956) 

450-590 


K 

336.75 

Macleod 
(1938) 

352-626 

Rohlin and 
Lodding (1962) 

350-490 

10. 

Li 

452 

Landolt- 
Bomstein 
tables (1969) 

465-605 

Murday and 
Cotts (1968) 

465-605 

li. 

Na 

371.05 

Macleod 

(1938) 

382-628 

i 

Meyer and 
Nachtrieb 
(1955) 

390-500 

■ 

Pb 


Landolt- 
Bomstein 
tables (1969) 

612-955 

Rothman and 
Hall (1956) 

630-810 

13. 

Rb 

'311 

Landolt- 
Bomstein 
tables (1969) 

375-465 

Norden and 
Lodding (1967) 

375-465 

14. 

m 


Walsdorfer 

(1983) 

530-660 

Careri et al 
(1958) 

530-660 


T1 

576 

Walsdorfer 

(1983) 


Cahill and 
Grosse (1965) 

525-640 

16. 

Zn 

692.6 

Yao and 
Kondic (1952) 

720-860 

Nachtrieb and 
Fraga (1963) 

700-860 
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Table 4.2 Hard sphere diameter parameters calculated from experimental viscosity data 


S. No 

Metal 

A 

BX10 4 

CX10 8 

Correlation 



(A 0 ) 

(A°/K) 

(A°/K 2 ) 

coefficient 

1 

Ag 

2.7471 

-1.5331 

1.0478 

0.999976 


Bi 

3.4504 

-8.2745 

2.9137 

0.99999 

3 

Cd 

2.8956 

-2.4961 

1.6283 

0.999995 

4 

Cs 

5.0116 

-7.4249 

-33.810 

0.99999 

5 

Cu 

2.0758 

4.182 

-20.164 

0.999984 

6 

Ga 

2.79 

-7.316 

33.262 

0.999999 

7 

Hg 

3.009 

-11.1 

73.418 

0.9998 

8 

In 

3.1111 

-5.735 

11.976 

0.99974 

9 

K 

4.4058 

-9.8522 

48.7957 

0.99984 

10 

Li 

2.8407 

-3.2767 

10.710 

0.999991 

11 

Na 

3.5624 

-6.21 

7.634 

0.999898 


Pb 

3.2811 

-4.0772 

6.1966 

0.999898 

13 

Rb 

4.7368 




14 



-2.4508 

-25.117 

0.999995 

15 

T1 

3.2257 

-4.496 

4.2375 

0.99998 


Zn 

2.65 

-2.9242 

11.19 

0.99998 
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4.3 Diffusion coefficients 

4.3.1 Hard sphere model 

The hard sphere diameters calculated by the above method were used to obtain the 
self-diffusion coefficients using the corrected Enskog expression. Furthermore, self- 
diffusion coefficients were also calculated using the Protopapas and Parlee (1974) hard 
sphere diameters. The resulting values are compared with the experimental values and 
these are plotted in figures 4.1 to 4.16. The computed values of self-diffusion coefficients 
were also listed in the tables (4.3-4.18) along with their experimental values. The hard 
sphere model gives comparable predictions for most metals. The obtained values are 
within 20-25% of the experimental values which themselves may entail experimental 
uncertainty of the order of 20%. In Fig 4.1-4.16, solid line indicates the Corrected Enskog 
equation using the diameter calculated from experimental viscosity data, dash-dot-dash 
line indicates the corrected Enskog values using the Protopapas and Parlee diameter, 
broken line and dot line indicate the Stokes-Einstein equation using the diameter 
calculated form experimental viscosity data and Protopapas and Parlee diameters 
respectively and circle line indicates experimental diffusivity values. 

The corrected Enskog hard sphere model underestimates the values of diffusion 
coefficients in most cases. In the case of Silver and Zinc, at low temperatures the 
Protopapas and Parlee diameter yields good values and as the temperature increases the 
extent of underestimation also increases. This model overestimates the values of 
diffusivities in the cases of Bi, In and Li. In the case of Bismuth, the Protopapas and 
Parlee diameter does very well compared to the other choices of diameter. In general all 
these models predict a weaker temperature dependence of the diffusivity than is observed 
experimentally. 

Within the context of the hard sphere model, the Protopapas-Parlee diameter, except 
for Gallium, does not display any marked superiority over the diameters obtained from 
viscosity data. This is despite the fact that the Protopapas- Parlee diameters were obtained 
from a statistical analysis of diffusivity data (Protopapas et al 1973). 



4.3.2 Stokes - Einstein Model 


An examination of tables (4.3-4.18) or plots (4.1-4.16) reveals that in most cases, the 
predicted values are within 5-25% of the experimental values. From these tables, it is 
readily seen that this model consistently underpredicts the value of diffusion coefficient, 
though the extent of underprediction varies one metal to another and/or with temperature. 
Thus for Ag, K, Li and Pb the predictions become progressively worse as the temperature 
increases whereas in the case of In, the predictions become close to the experimental 
values. Yet in some other systems e.g. Sn and Rb, the divergence between the 
experimental and predicted values decreases as the temperature increases. Suffice it to 
say here that similar deviations have been widely reported in the literature even for 
ordinary liquids. 

Some investigators have attributed these deviations to the non-spherical shape of 
diffusing molecule, and/or to solute-solvent interactions, etc., while numerous 
investigators have suggested ad hoc empirical modifications to improve the 
correspondence between the predictions of this model and experimental values of self- 
diffusion coefficient. Based on extensive data and heuristic considerations, many workers 
(Tyrrell, 1981; Pollack, 1981; Pollack and Enyeart, 1985) have argued that the self- 
diffusion coefficient depends upon the medium viscosity as: 

Doc-L (4.2) 

P 

where scores of values of p ranging from 0.53 to 0.69 have been documented in the 
literature and indeed Ziwanzig and Harisson (1985) have argued that there are no 
theoretical reasons to expect a single value of p to be universally applicable to all 
systems. Chauhan e t al (2000) have attempted this modification in the study of liquid 
metals, and the values of p ranging from 0.62 to 1.56 were obtained, showing absolutely 
no correlation with any structural, molecular or thermodynamic quantity. This more or 
less reinforces the notation of Ziwanzig and Harisson (1985). It is worthwhile to point out 
here that such a wide variation in the value of p cannot be attributed to experimental 
uncertainty in the diffusion coefficient or viscosity. The only other way to improve the 
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Table 4.3 Self-diffusion coefficients and hard sphere diameters for liquid Ag 


Temp 

(K) 

^Calculated 

(A 0 ) 

GProtopapas 

(A 0 ) 

D£f‘xl0>V S ) 

D S _ E xl0 9 (m 2 /s) 

A&/»xl0 9 (m 2 /s) 

^Calculated 

Gprotopapas 

^Calculated 

Gprotopapas 

1260 

2.5706 

2.5828 

3.1093 

2.7962 

2.8872 

2.8735 

2.7214 

1280 

2.5680 

2.5802 

3.2127 

2.8962 

3.0270 

3.0127 

. 

2.8574 

U 







2.9901 

a 


2.5750 





3.1275 

1340 

2.5605 

2.5725 

3.5318 

3.2037 

3.4486 

3.4326 

3.2668 



2.5699 

3.6411 




3.4080 






3.7428 

3.7254 

3.5509 


2.5530 

2.5650 

3.8638 

3.5220 

3.8937 

3.8756 

3.6950 


Table 4.4 Self-diffusion coefficients and hard sphere diameters for liquid Bi 


V 

1 


Oprotopapas 

(A 0 ) 

D En S k° g xl0 9 (m 2 /s) 

D S _ E xl0 9 (m 2 /s) 


^Calculated 

OprotoDauas 

^Calculated 


540 

3.0885 

3.1466 

2.2484 

1.5239 

2.0926 

2.0539 

0.8002 

560 


3.1394 

2.4347 

1.6448 

2.2915 

2.2470 

0.8700 

580 

3.0685 

3.1322 

2.6246 

1.7687 

2.4973 

2.4465 

0.9403 

m 

3.0588 

3.1252 

2.8178 

1.8956 

2.7094 

2.6518 

1.0111 

620 

3.0494 



2.0254 

2.9275 

2.8627 

1.0822 

640 

3.0402 

3.1116 

3.2123 

2.1580 

3.1517 

3.0794 

1.1533 


3.0312 

3.1049 



3.3813 

3.3011 

1.2244 


3.0225 

3.0983 

3.6147 

14311 

3.6166 

3.5280 

1.2953 
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Table 4.5 Self-diffusion coefficients and hard sphere diameters for liquid Cd 


Temp 
( K) 

^Calculated 

(A 0 ) 

^Protopapas 

(A 0 ) 

DS?x 10 5 (m ! /s) 

D s _ e xl0 9 (m 2 /s) 

7)&pixl0 9 (m 2 /s) 

^Calculated 

^ProtODauas 

^Calculated 

CT Protopapas 

610 






1.9642 

2.3595 

m 

2.7448 

2.7484 

2.2559 

2.1916 

2.1361 

2.1333 

2.5732 

a 




2.3391 

2.3096 

2.3075 

2.7913 

670 

2.7357 

2.7372 

2.5191 

2.4899 

2.4884 

2.4870 

3.0131 

■ 



2.6557 

2.6438 

2.6719 

2.6713 

3.2383 

u 




2.8007 



3.4462 

730 

2.7221 

2.7210 

2.9389 


3.0527 


3.6963 


2.7176 

2.7158 



3.2496 

3.2518 

3.9282 


Table 4.6 Self-diffusion coefficients and hard sphere diameters for liquid Cs 


I 

^Calculated 

(A°) 

OProtopapas 

(A 0 ) 

D£* xl0 9 (m 2 /s) 

D s _ e xl0 9 (m 2 /s) 

DsxptX 10 9 (m 2 /s) 

^Calculated 

O’ProtODapas 

^Calculated 

OProtopapas 

325 

4.7346 

4.7398 

2.6665 

2.6045 

2.5258 

2.5230 

2.6945 

335 

4.7249 


2.8520 

2.7864 

2.7488 

2.7458 

2.9364 

fi| 

4.7152 

4.7208 

3.0446 

2.9727 

2.9781 

2.9746 

3.1814 

355 







3.4370 

365 

4.6955 



3.3577 

3.4602 

3.4552 

3.6945 

375 



3.6651 

3.5563 

3.7131 

3.7069 

3.9559 

385 

4.6756 

4.6846 

3.8863 

3.7586 

3.9721 

3.9645 

4.2209 

395 

4.6656 

4.6758 

4.1145 1 

3.9647 

4.2386 

4.2293 , 

4.4888 
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Table 4.7 Self-diffusion coefficients and hard sphere diameters for liquid Cu 


Temp 

(K) 

^Calculated 

(A 0 ) 

Gprotopapas 

(A 0 ) 

Dg* xlO s (m ! /s) 

D s _ e xl0 9 (m 2 /s) 

-D&^xl0 9 (m 2 /s)' 

^Calculated 

^Protooaoas 

^Calculated 

G’ProtopaDas 





3.2933 

3.8456 

3.8283 

4.1231 

KH 



3.7170 

3.3492 

3.9192 

3.9011 

4.2311 



2.2784 

3.7889 

3.4055 




m 



3.8629 

3.4622 




1410 

2.2646 

2.2763 

3.9391 




4.5622 

1420 

2.2631 

2.2753 

4.0176 



4.1976 

4.6749 

1430 







4.7888 



2.2732 

4.1812 

3.6922 

4.3746 

4.3489 

4.9038 


Table 4.8 Self-diffusion coefficients and hard sphere diameters for liquid Ga 


Temp 

(K) 

^Calculated 

(A 0 ) 

Gprotopapas 

(A 0 ) 

DS? X 

I0 9 (m 2 /s) 

D s _ e xl0 9 (m 2 /s) 

f5 £ ^xl0 9 (m 2 /s) 




^Calculated 

O’Protopauas 

^Calculated 

CProtopauas 


330 

2.5848 

2.5622 

1.5925 

1.9541 

1.4417 

1.4544 

1.9937 

335 

2.5822 

2.5596 



1.4990 

1.5123 

2.0452 

l^jjj 

2.5797 



2.0757 

1.5573 

1.5711 

2.0964 



2.5546 

1.7580 

2.1374 

1.6165 

1.6308 

2.1473 

a 

2.5747 

2.5521 

1.8145 

2.1997 

1.6764 

1.6913 

2.1981 

355 

2.5722 

2.5496 

1.8716 

2.2625 

1.7372 

1.7526 

2.2485 

360 

2.5697 



2.3259 

1.7988 

1.8148 

2.2986 

365 

2.5673 

2.5446 

1.9877 

2.3899 

1.8612 

1.8778 

2.3484 












































































Table 4.9 Self-diffusion coefficients and hard sphere diameters for liquid Hg 


Temp 

(K) 

^Calculated 

(A 0 ) 

CProtopapas 

(A 0 ) 

10> ! /s) 

D s _ e xl0 9 (m 2 /s) 

D Exp ix 10 9 (m 2 /s) 

^Calculated 

OProtoeapas 

^Calculated 

^Protopapas 

270 

2.7628 

2.7722 

1.3106 

1.2250 

1.2626 

1.2583 

1.2938 

■ 

2.7488 

2.7584 



1.4722 


1.4732 

310 

2.7355 

2.7450 




1.6865 

1.6496 

330 

2.7227 



1.7848 


1.9147 

1.8219 

350 

2.7104 

2.7196 

2.0905 


2.1579 

2.1506 

1.9890 

370 

2.6988 

2.7075 

2.2939 

2.1891 

2.4026 

2.3949 

2.1520 


Table 4. 10 Self-diffusion coefficients and hard sphere diameters for liquid Pb 


Temp 
( K) 

^Calculated 

(A 0 ) 

CFprotopapas 

(A 0 ) 

Dg* xl0 9 (m 2 /s) 

D,. £ xl0 9 (m 2 /s) 

Dgxptx 10 9 (m 2 /s) 

O' Calculated 

0>Protopapas 

^Calculated 

Oprotopapas 

630 

3.0488 

3.0647 

1.9773 

1.7671 

1.7516 

1.7426 

2.6160 

U 

3.0423 


2.1059 

1.8847 

1.9146 

1.9045 

2.9189 

670 

3.0357 

3.0522 

2.2379 

■ 

2.0047 

2.0838 

2.0726 

3.2350 

690 

3.0293 

3.0462 



2.2594 

2.2469 

3.5640 

710 

3.0229 

3.0402 

2.5113 

2.2520 

2.4409 

2.4270 

3.9053 

730 

3.0165 

3.0343 

2.6526 

2.3791 

2.6282 

2.6128 

4.2578 




2.7968 

2.5083 

2.8211 


4.6207 

770 

3.0039 

3.0227 

2.9439 

2.6397 

3.0192 

3.0004 

4.9930 

790 

2.9977 

3.0171 

3.0936 

2.7731 

3.2228 

3.2021 

5.3749 

810 

2.9915 

3.0115 

m 


3.4312 

3.4085 

5.7645 
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Table 4. 12 Self-diffusion coefficients and hard sphere diameters for liquid Li 



525 

2.6982 

2.6898 

10.5965 

11.1997 

8.5676 

8.5942 

9.1969 

545 

2.6939 

2.6829 

11.7730 

12.6085 

9.2840 

9.3223 

10.1740 

565 

2.6898 

2.6760 




10.0699 

11.1740 

585 

2.6857 

2.6693 

14.2953 

15.6322 

10.7678 

10.8338 

12.1950 

605 


2.6627 

15.6383 

17.2423 

11.5315 

11.6136 

13.2320 








































































Table 4.13 Self-diffusion coefficients and hard sphere diameters for liquid Na 


Temp 

(K) 

^Calculated 

(A 0 ) 

GProtopapas 

(A 0 ) 

D£f‘xl0> ! /s) 

D S _ E xl0 9 (m 2 /s) 

£>Ex P tX 10 9 (m 2 /s) 

^Calculated 

OProtODauas 

^Calculated 

^Protopapas 

390 

3.3313 

3.3246 

4.4786 

4.6857 


4.0226 

4.7828 

m 

3.3257 

3.3191 





5.1728 

E9 

3.3145 

3.3083 

5.4895 




5.9834 

mm 

3.3033 

3.2978 

6.2173 

6.4143 

5.7318 

5.7414 

6.8300 

460 

3.2922 

3.2875 



6.4929 

6.5022 

7.7070 

480 

3.2812 

3.2774 

7.7973 

7.9445 

7.2956 

7.3040 

8.6100 

500 

3.2702 

3.2676 

8.6467 

8.7545 

8.1359 

8.1425 

9.5330 


Table 4.14 Self-diffusion coefficients and hard sphere diameters for liquid In 


Temp 

(K) 

^Calculated 

(A 0 ) 

Oprotopapas 

(A 0 ) 

10 9 (m 2 /s) 

D s _ e xl0 9 (m 2 /s) 

T»£xp/xl0 9 (m 2 /s) 

O' Calculated 

OfProtopapas 

^Calculated 

Cprotopapas 

450 

2.8773 

2.8935 

2.1291 

1.8847 

1.9035 

1.8928 

1.8890 

470 


2.8853 

2.3340 

2.0620 

2.1213 

2.1087 

2.1225 


2.8588 

2.8772 


2.2446 

2.3478 

2.3328 

2.3612 


2.8498 

2.8693 



2.5827 

2.5651 

2.6046 

i 

530 

2.8408 

2.8616 

2.9936 

2.6244 

2.8254 

2.8049 

2.8526 




3.2273 



3.0519 

3.1032 

570 

2.8231 

2.8465 

3.4673 

3.0225 

3.3332 

3.3058 

3.3559 




3.7134 



3.5663 

■ 

r 

3.6100 






























































Table 4.15 Self-diffusion coefficients and hard sphere diameters for liquid Rb 


Temp 

(K) 

^Calculated 

(A 0 ) 

Gprotopapas 

(A 0 ) 

D£!°* xlO s (m’/s) 

D s _ e xl0 9 (m 2 /s) 

Z)^,xl0 9 (m 2 /s) 

^Calculated 

CProtODanas 

^Calculated 

GProtopapas 

375 

4.3643 

4.3672 

3.9896 

3.9410 

3.9259 

3.9234 

4.6304 

385 

4.3554 

4.3591 

4.2401 

4.1764 

4.2151 

4.2116 

4.9612 

395 

4.3466 

4.3511 

4.4970 

4.4166 

4.5136 

4.5090 

5.2972 

405 

4.3378 

4.3432 

4.7600 

4.6614 

4.8193 

4.8133 

5.6376 



4.3354 

5.0291 

4.9107 

5.1340 

5.1266 

5.9819 

425 

4.3204 

4.3277 

5.3040 

5.1643 

5.4559 

5.4468 

6.3296 

445 

4.3033 

4.3125 

5.8706 

5.6840 

6.1219 

6.1087 

7.0329 

465 

4.2863 

4.2977 

6.4586 

6.2195 

6.8174 

6.7993 

7.7439 


Table 4,16 Self-diffusion coefficients and hard sphere diameters for liquid Sn 


Temp 

(K) 

^Calculated 

(A 0 ) 

^Protopapas 

(A 0 ) 

Dg?*x 10 9 (m 2 /s ) 

D s _ e xl0 9 (m 2 /s) 

DExpt* 10 9 (m 2 /s) 

^Calculated 

Ctyrotopapas 

^Calculated 

0>Protopapas 

530 

2.9288 

2.9313 

2.0913 

2.0506 

1.9253 

1.9236 

2.6070 



2.9242 




2.1580 

2.8498 

570 

2.9079 

2.9172 

2.5455 

2.3846 



3.0900 



2.9104 

2.7943 

2.5578 

2.6757 

2.6636 

3.3440 

610 

2.8862 

2.9037 

3.0579 

2.7350 

2.9519 

2.9342 

3.5947 

630 

2.8751 

2.8970 

3.3365 

2.9159 

3.2411 

3.2165 

3.8460 


2.8638 

2.8905 

3.6303 

3.1005 

3.5423 

3.5096 

! 4.0970 

660 

2.8580 

2.8873 

3.7829 

3.1941 

3.6976 

3.6601 

4.2230 


41 





























































Table 4.17 Self-diffusion coefficients and hard sphere diameters for liquid T1 


Temp 

(K) 

^Calculated 

(A 0 ) 

OProtopapas 

(A 0 ) 

Dgr* io> 2 /s) 

D S _ E xl0 9 (m 2 /s) 

D Expt x\(P(vci'Js) 

^Calculated 

^ProtODaoas 

^Calculated 


525 

3.0013 

3.0146 

1.4081 

1.2592 

1.2140 

1.2086 

1.3925 

550 

2.9912 

3.0061 

1.5681 

1.3921 

1.3950 

1.3881 

1.6868 

600 

2.9712 

2.9897 

1.9124 

1.6715 

1.7904 

1.7794 

2.3591 

610 

2.9672 

2.9865 

1.9850 

1.7295 

1.8747 

1.8626 

2.5063 

620 

2.9632 

2.9833 


1.7881 



2.6574 

630 

2.9593 

2.9801 

2.1338 

1.8474 



2.8124 

640 

2.9553 

2.9770 

2.2101 

1.9073 

2.1373 

2.1217 

2.9711 


Table 4. 1 8 Self-diffusion coefficients and hard sphere diameters for liquid Zn 


Temp 

(K) 

^Calculated 

(A 0 ) 

0*Protopapas 

(A 0 ) 

D“" f xl0’(m 2 /s) 

D s _ e xl0 9 (m 2 /s) 

A&pfxl0 9 (m 2 /s) 

^Calculated 

GProtopapas 

^Calculated 

GProtouapas 



2.4572 





2.1411 



2.4528 

1.6435 

2.6581 


1.4973 

2.3944 


2.4949 

2.4484 

1.7334 

2.8201 

1.5943 

1.6246 

2.6616 

780 

2.4900 

2.4399 

1.9147 

3.1540 

1.8552 

1.8933 

3.2357 

800 

2.4877 

2.4357 



1.9918 

2.0344 

3.5416 

820 

2.4855 

2.4315 



2.1325 

2.1798 

3.8594 

840 

2.4833 

2.4274 

2.1884 

3.6782 


2.3295 

4.1885 

860 

2.4813 

2.4234 

2.2795 

3.8588 



4.5284 
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Diffusion coefficient, 10 m/s 



Temperature, K 

Fig 4.3 Self-diffusion coefficients of liquid Cadmium as a function of temperature. 



Temperature, K 


Fig 4.4 Self-diffusion coefficients of liquid Caesium as a function of temperature. 
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Diffusion coefficient, 10" 9 m 2 /s 


1360 1370 1380 1390 1400 1410 1420 1430 1440 1450 1460 

Temperature, K 

Fig 4.5 Self-diffusion coefficients of liquid Copper as a function of temperature. 



Temperature, K 


Fig 4.6 Self-diffusion coefficients of liquid Gallium as a function of temperature^ 
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Diffusion coefficient, 10" 9 m 2 /s 



Fig 4.8 Self-diffusion coefficients of liquid Indium as a function of temperature. 
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440 460 480 500 520 540 560 580 600 620 

Temperature, K 

Fig4,10 Self-diffusion coefficients of liquid Lithium as a function of temperature. 





Diffusion coefficients, 10' 9 m 2 /s 



380 400 420 440 460 480 500 520 

Temperature, K 

Fig 4.1 1 Self-diffusion coefficients of liquid Sodium as a function of temperature. 



620 640 660 680 700 720 740 760 780 800 820 

Temperature, K 

Fig 4, 12 Self-diffusion coefficients of liquid Lead as a function of temperature. 
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520 540 560 580 600 620 640 660 680 

Temperature, K 

Fig 4. 14 Self-diffusion coefficients of liquid Tin as a function of temperature. 




X 


> I i ' 1 1 1 i 1 1 T T i » • 1 * I * • 

680 700 720 740 760 780 800 820 840 860 880 

Temperature, K 

Fig 4. 16 Self-diffusion coefficients of liquid Zinc as a function of temperature. 





predictability of this model is to reconsider the assignment of a suitable value to size of 
the diffusing species, R. So, in this regard, Goldschmidt atomic diameter was replaced by 
the effective hard sphere diameters, which were calculated using the approach outlined 
by Protopapas and Parlee and from the experimental viscosity data. 

It was observed that the predictions of the diffusion values by these two diameters 
are almost identical. This points to the reduced sensitivity the Stokes-Einstein model to 
the choice of the characteristic diameter. It was shown earlier (Chauhan, 2000) that 
within the context of the Stokes- Einstein model, the temperature-independent. Gold 
schmidt diameter consistently under predicts the diffusion coefficient and the use of the 
temperature dependent Protopapas-Parlee diameter and diameter calculated from the 
viscosity data improves predictions. Thus, while incorporating the temperature 
dependence of the characteristic diameter helps in improving the Stokes-Einstein model 
predictions, neither of the two diameters considered here can capture the temperature 
dependence of diffusivities adequately. 

Finally, it is appropriate to point out here that for estimating the Protopapas and 
Parlee diameter, one needs only the melting point temperature and the corresponding 
value of density. On the other hand, if one were to use the viscosity data as illustrated in 
this work, the parameters A, B and C have to be determined for each metal. However, it 
is worth nothing that the Protopapas and Parlee diameter was obtained by using 
diffusivity data and hence is to be regarded more as a fit to the diffusivity data rather than 
as a method for a priori determination of diffusivity of liquid metals. 

4.4 Corresponding states model 

Chapman (1966) has shown that it was possible to collapse viscosity data onto a 
single curve by the adjustment of one microscopic parameter (e). It is found that this 
empirically evaluated parameter has some fundamental significance. Pasternak (1967) 
has obtained a correlation for self-diffusion coefficient and viscosity of liquid metals 

tRlo jo TMTT 

ns 

w»ifc,-'-.3To B\ 1 3 3 B R f) 
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Fig 4.17 


Plot of fj.*{v*y vs. T' using long-range oscillatory potentials 
(Johnson et al., 1964). 



Fig 4.18 Plot of r 


(Johnson et al., 1964). 


vs. T* using long-range oscillatory potentials 
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Fig 4.19 Plot of vs. (r*y using long-range oscillatory potentials 
(Johnson et al., 1964). 
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Fig 4.20 Plot of —■ vs. V using long-range oscillatory potentials 
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(Johnson et al., 1964) 
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Fig 4.25 Plot of i u* (f* ) 3 vs. T* using Lennard- Jones (m-n) potential parameters 
(Zhen and Davies, 1983). 



D 

Fig 4.26 Plot of - vs. V using Lennard-Jones (m-n) potential parameters 

(V) 

(Zhen and Davies, 1983). 
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using Helfand and Rice (1960) model and the parameters, which were obtained by 
Chapman (1966). The Chapman’s correlation was tested with various potential 
parameters, which were obtained from different types of potentials. Lennard-Jones (m-ri) 
(Zhen and Davies, 1983), long-range oscillatory (Johnson et al., 1964) and pseudo 
potentials (Ranganathan and Pathak, 1994) were chosen in this approach. Different 
combinations of reduced quantities were calculated and were plotted against reduced 
temperatures. The trends of such plots are shown in figures (4.17-4.28). 


We restricted our selves to certain metals for which both viscosity and diffusivity 
data is available. The plots of volume dependency of reduced viscosity and reduced 
diffusivity (fig 4.17 and 4.18) using long-range oscillatory potential parameters, shows 
that liquid metals do fall on the single curve, but the plots of ratio of reduced viscosity 
and reduced diffusivity (fig 4.19-4.20) show that the liquid metals Hg, Cs, Pb do fall onto 
one curve where as K and Na follow another curve. The reduced quantities, calculated 
using the pseudo potential parameters (Ranganathan and Pathak, 1994) were shown (fig 
4.21-4.24) to follow the same trend irrespective of combinations of reduced quantities. 
The curves obtained using pseudo potentials are close to each other than those obtained 

/ \1 £)* 

from the other potentials. The plots of jU*(F* j 3 vs. T* and — vs. T* (fig 4.25 and 

(v} 

4.26) show that the liquid metals are not falling on the same single curve. But the plot of 
vs. (t”)* shows that the liquid metals Zn, Na, and K collapse onto the same line 


where as liquid Rb tends to fall on that curve. The plot of 




vs. J* shows that only 


Na and K are falling on the same line and the rest are falling on different curves. 


All these potentials were used in the calculations of bulk modulus, electrical 
resistivity and the other thermodynamic properties of liquid metals. These parameters 
lead to the values which show considerable agreement with experimental values. 
However, it is not possible to pinpoint the exact reason for this type of behavior. 
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Furthermore, no potential gives potential parameters for all liquid metals. The 
temperature range, which we used in our calculations, is +10% of the melting 
temperature. So this temperature range changes from one metal to another, therefore, the 
reduced quantities range also varies from one metal to the other. This might be the reason 
why these metals do not follow the same curve, because according to definition of 
corresponding state principle, the reduced viscosity or diffusivity of similar substances in 
any thermodynamic state should be a universal function of the reduced thermodynamic 
properties of that state such as reduced temperature and pressure. One way to obtain such 
universal function is to follow the scheme outlined by Chapman (1966) and Pasternak 
(1967) with different types of potential parameters and then minimize the errors in 
correlating a universal function. 
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Chapter 5 


CONCLUSIONS 


5.1 Conclusions 


In this work, the efficacy of a hard sphere diameter, inferred from viscosity- 
temperature data, in predicting the self-diffusion coefficient of liquid metals using the 
corrected Enskog and the Stokes-Einstein relationship has been extensively investigated. 
An attempt has been made to examine the applicability of the corresponding state 
principles for a universal correlation for viscosity and diffusivity of the metallic system. 
Based on the analysis reported herein, following conclusions can be drawn: 

1. It has been shown that hard sphere diameter estimated from viscosity data can lead to 
reasonably accurate values of self-diffusion coefficients of liquid metals. The relative 
ease with which such viscosity data is available makes this approach attractive. However, 
with the above hard sphere diameters, neither the Stokes-Einstein equation nor the 
corrected Enskog equation are able to predict the temperature dependence of the 
diffusivity as strongly as observed experimentally. Nevertheless, the predictions are 
comparable to those obtained using the Protopapas and Parlee diameters which are 
inferred from diffusion data itself. 

2. Liquid metals do not follow the universal curve of corresponding states model when 
tested with various types of potential parameters. The main difficulty appears to be in the 
lack of information regarding the potential parameters for the systems investigated 
herein. 


5.2 Recommendations for future work 


Further experimental work is needed to encompass many more metallic systems in 
terms of their viscosity and diffusivity data over a wide range of temperature as possible. 
This will enable the validation of the approaches employed in this work. Similarly, in 
order to test the applicability of the corresponding state principle in a much more through 
manner, it is required to evaluate the potential parameters for molten metals so that a 
coherent framework can be developed. 
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APPENDIX A 



Fig A: Flow chart for the calculation of hard sphere diameter from 
experimental viscosity data. 
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APPENDIX B 


Table B Liquid metals and their physical properties. 


S.No 

Metal 

M 

T m 

(K) 

Pm (10 3 'kg/m 3 ) 

-A (10' 1 kg/m 3 K) 

1 

Ag 

107.87 

1233.65 

9.33 

9.10 i 

2 

Bi 


544.1 

10.05 

11.8 : 

3 

Cd 

112.24 

594 

8.01 

12.2 , 

4 

Cs 

132.9 

301.65 

1.84 

5.70 j 

. 

5 

Cu 

63.54 

1356 

8.00 

8.01 ; 

6 

Ga 

69.72 

303 

6.10 

5.60 | 

' 

7 

Hg 

200.5 

234.28 

13.69 

24.2 } 

' 

■ ' I 

8 

In 

114.82 

429.55 

7.03 

6.80 j 

!■ 

9 

K 

39.09 

336.75 

0.83 

2.40 j 

' ' ' 

10 

Li 

6.94 

452 

0.518 

1.00 

11 

Na 

22.98 

371.05 

0.927 

2.35 ! 

12 

Pb 

207.2 

600.55 

10.67 

13.2 

13 

Rb 

85.5 

311 

1.48 

4.50 

mm 

Sn 

118.6 

505 

6.98 

6.10 


T1 

204.3 

576 

11.35 

13.0 

E 

Zn 

65.38 

692.6 

6.58 

9.80 















































